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I. INTRODUCTION 



A method is developed whereby without any undue ne- 
glect of certain factors the decrease in the mass flow of 
air and in the drag of the heated radiator, i.e., one un- 
der actual operating conditions, as compared with the cor- 
responding values in the cold state, may be simply computed. 
Although a symmetrical duct radiator has been used to 
bring out the flow relations, the results apply equally as 
well to the unsymmetr ical radiator shapes usual in air- 
plane construction. 

It was ivieredith (reference l) who first pointed out 
the possible gain in power required to tow the radiator 
system under the operating condition as compared with the 
cold condition and showed how the momentum or power gain 
resulting from the heating of the air by the coolant may 
be approximately computed with the aid of simple momentum 
and energy considerations. Since the computation, however, 
involves rather extensive simplifications, the results may 
be used tentatively as a first approximation only. 

A very clear presentation of the effect of the air 
forces on a duct radiator, particularly on the radiator 
cowling, has been given by Weise (reference 2) who ex- 
plained the change of state of the air in the radiator 
with the aid of a velocity-pressure diagram. In the pres- 
ent paper simple formulas are derived with the aid of 
which the lowering in the drag of a radiator due to heat- 
ing may be computed with an accuracy sufficient for prac- 
tical requirements. There is first considered the case 
of the cold radiator and next the effect of the heating 
of the radiator on the rate of air flow and on the drag. 



♦"Beitrag zur Theorie des beheizten Dusenkuhlers . " From 
Luftfahrtf orschung, vol. 15, no. 10/n, Oct. 10, 1938, pp. 
500 to 504. 
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II. NOTATION 



v o- p o- F 0' P b » velocity, air density, cross-sec- 
tional area, and pressure at 
duct inlet. 

v i» ^ • Pi» corresponding values immediately 

in front of radiator "block. L See 

/fig. 1 

v 2 , P 3 , F, p 2 , corresponding values "behind radia- 
tor "b lock. 



v , P , F , p QI corresponding values at duct exit 

p = ~. p; * p 0 . p 3 * p 2 

Wj-jj., internal drag of cold radiator. 

W i , internal drag of heated radiator. 

c w , drag coefficient of cold radiator. 

cpjj., effective rate of air flow v x /v 0 (equati on (6) ) 
through cold radiator. 

cp 0 , theoretical flow ratio through cold radiator (equa- 
tion (6a) ) . 

cp flow ratio through the heated radiator (equations 

(14) and (17) ) . 

T|j-) , efficiency of diverging passage ahead of radiator 
. 1 "block (equations (7) and (7a)). 

v 3 

CD„ = . • 

° ' V ° 

A cp = 1 i. 

cPk , 

N, engine power (hp.). 

Q = 75 a N, heat quantity to he conducted away, mkg/s. 



N.A.C.A. Technical Memorandum . No . 893 



3 



6 = 



75 ccN 



— P V 3 F 

2 H o v o * 



< n t k« thermodynamic efficiency of the radiator elements 

6, temperature drop between cooling fluid and air 
radiator" entry. 



K = 1.4; 0-0 = 0. 24; A = — L_ (kcal/kgm). 
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III. THE UNHEATED DUCT RADIATOR 



We first consider the flow through and around an un- 
heated duct radiator. On figure 1 the cross sections and 
the corresponding velocities, densities, and pressures 
are indicated. For the cold radiator Vj = v a . The cross- 
sectional area F 0 , which is proportional to the quantity 
of air flowing through, depends as is known on the cross- 
section ratio E/F a , on the drag of the enclosed radiator 
blo.ck,. and on the remaining flow losses within the radia- 
tor.. Of these losses that, due to the widening of the 
flow cross section from E 0 to F, is of particular im- 
portance. The flow may generally be considered only part- 
ly as a flow through a diffuser in the usual sense, i.e., 
as a flow through a diverging passage, since the flow 
cross section' begins to enlarge ahead of the duct entry 
itself. Up to the stagnation point at the nose of the 
radiator, the flow may be considered as being free from 
losses but in the diverging duct in front of the radia- 
tor elements losses arise and may be taken approximately 
into account by introducing a diffuser efficiency as will 
be shown below. However, losses are also involved in the 
external flow about the radiator. This flow may be com- 
pared with "that about the upper side of a more or less 
steeply inclined wing section. The less the air quantity 
flowing through for given dimensions of the duct radiator 
and given speed, i.e., the more divergent the stream lines 
ahead of the radiator, the greater are the losses of the 
external flow, particularly those about the leading edge 



4 



N.A.C.A. Technical Memorandum. No. 893 



of the radiator. The flow pattern then corresponds ap- 
proximately to that about ah airfoil set at a large angle 
of attack; i.e., the flow may, under certain conditions, 
separate from the outside duct wall and give rise to an 
external radiator drag that may amount to a considerable 
portion of the total drag. 

We shall next consider the flow at the radiator exit 
section. Since the flow at the outer duct wall does not 
entirely separate and the air exit from the radiator takes 
place without appreciable contraction, the static pressure 
at the plane of the exit section may be assumed as equal 
to the external pressure p Q . Not considering for the mo- 
ment the external flow losses, there is obtained the veloc- 
ity distribution shown by the continuous line in figure 1. 
From the momentum theorem, neglecting changes in density, 
the "internal drag" of the cold radiator may be computed as 



W, 



■k 



= P„ F v 



(v 0 - v 3 ) 



or 



W 



= P 0 F v 0 < 



cp k (1 - p cp k ) 



(1) 



where cp^ = — and (3 = 



F 



The internal drag thus depends on the opening ratio 
of the exit passage and on the flow ratio cp^. The latter 
may be readily computed if the drag coefficient (c w ) of 
the cold radiator elements is known, for -example from 
tests in the closed duct. The pressure rise ahead of the 
radiator elements not considering diffuser loss is ob- 
tained from the Bernoulli equation 



- P, 



_ P o / ; 
- — (v 0 



- V 



(2) 



For simplicity the diffuser loss is expressed in terms of 
the. total theoretical pressure rise although, as previously 
mentioned, the diverging flow only partly takes place 
through the duct itself, and for the actual pressure rise 
there is obtained 



- v 



) 



\ 



(3) 



T. 
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" The- f low- in" the duct 'behind the radiator may be con- 
sider ed: as pract.icall-y ..without losses, so that 

Pa = Po + ~ (j3 a , 1) (4) 

For the pressure drop through the radiator elements, we 
have 

Po o 

Pi - Ps = -jT v i c w (5) 



and after simple substitutions, there is obtained from 
equations (3), (4), and (5) the rate of flow ratio for the 
cold radiator 

= ~ = / ^r—z — — - (e) 

v o v c w + p - (1 - r\ Di ) 

This ratio thus depends on the diffuser efficiency as well 
as on c w and. j3 . for .otherwise equal conditions this 
efficiency will be .somewhat lower for a duct radiator with 
small, inlet section than for one with large inlet section. 
A contraction of the inlet opening thus leads to a lower- 
ing in' the air discharge rate (equation (6)) but generally 
to an increase in the internal drag (equation (l)). 

The actual velocity distribution at the radiator exit 
corresponds approximately to the dotted curve in figure 1. 
The additional momentum loss is due to the external radia- 
tor loss, i.e., to the frictional resistance of the exter- 
nal cowling. This loss will be relatively small, provided 
the flow separation is not too great, but will increase 
with decreasing discharge< ratio (i.e., for large c w and 

on account of the strong flow about the leading edge 
and, since the internal, drag is relatively- small, it will 
form a considerable portion of the total drag. For the 
same reason, this loss will be somewhat smaller for a ra- 
diator with diffuser than for one without a diffuser. The 
inlet diffuser thus leads to a lowering of the external 
drag at the expense of the internal drag and the rate of 
discharge through the radiator. These changes will, how- 
ever, not be very large, so. that, the total drag may be 
considered as sufficiently independent of the inlet cross 
section while it depends very much on the value c w of 
the. enclosed radiator block and on the ratio,- radiator 
cross section/exit cross section. 



6 



N.A.O.A.. Technical Memorandum -No . 893 



If the rate of flow ratio of the cold radiator is 
known from measurements, the diff user : eff i ci ency may con- 
versely be determined from it. Denoting by 

(6a) 




the theoretical discharge ratio of the duct radiator is 
obtained from equations (6) and (6a): 

/ ca_ 3 1 - cp> 3 



or, since cp Q and cp^ are, in general, not very differ- 

= (sr) < 7a > 



ent 

^D 



a very simple relation for the expansion efficiency. This 
efficiency may be de t ermined , f or example, from the effec- 
tive rate of discharge ratio cp^. obtained from model 

tests and from the theoretical discharge ratio <P 0 , 



IV. THE HEATED DUCT RADIATOR 

The heating of the radiator elements leads to a de- 
crease in the air density from P x to P 2 and to an in- 
crease in the velocity from v x to v g . (See fig. 1.) 

Eor the two-control surfaces aa and bb directly 
in front of and behind the radiator block, the momentum 
equation gives 

Px + P o v i 2 = Pa + p 2 v 2 2 + f < 8 )- 

where W denotes the total frontal and frictional resist- 
ance of the radiator block 1 . For the cold radiator (v x = 

1 W is, of course, not the equvalent to the internal drag 
of the system: cooler block + duct, according to equation 
(1). 
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v s ) , th,e pressure diff erence is a measure, of this resist- 
ance: 

, W = F.( Pi - p 2 ) = Z° Cw ? 

C 2 • ; 

For the heated; radiator., this is no . longer the case since, 
as follows from equation (8), the change in momentum here 
also plays an important part. For equal weight of cooling 
air, the pressure drop in the heated radiator is greater 
than in the. cold one since the momentum at the exit ex- 
ceeds that at the inlet and moreover the mean dynamic pres- 
sure increases to 



In the drag coefficient c 



w • 



there is included the 
frontal resistance, the frictional loss, for example, in 
the cooling elements, and the mixing loss at the radiator 
exit. The frictional drag depends on the Bey.noldB Number, 
The small difference in c w due to the difference in the 
discharge rate and in the mean air temperature obtained 
with the heated radiator as compared with the cold radia- 
tor, will not, however, be taken into account in the fol- 
lowing computation, since, as further computation shows, 
the change in the discharge rate is very slight and the 
frictional resistance constitutes only a (not accurately 
known generally) portion of the total drag. 

On substituting the mean dynamic pressure, (equation 
(8) goes over into 3 



Ap h = Pj - p a = (p 3 v 2 2 



- Pov.a) +'|^ Vi a + ^a*) Cw (8a) 



If the drag coefficient of the hot radiator is defined as 

r ^h 
Cw h - n 



Ap-^ includes the pressure rise due to the momentum in- 
crease and the inoref.se in the mean dynamic pressure. The 
drag coefficient will then of course differ essentially 
(Continued on next page) 
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On account of the increased pressure drop in the radiator 
"block, the streamlines ahead of the radiator must diverge 
more strongly than in the ease of the cold radiator; i.e., 
there is a decrease in therate of air discharge until an 
equilibrium condition corresponding to equation (8a) is 
reached. We shall therefore first determine the cooling 
air discharge rate through the heated radiator and from it 
the internal drag of the heated radiator. 

Equation (8a) contains five unknowns (p x , p 2 , v 1 , v 2 
p 3 ) . If we neglect the very small change in the diffuser 
efficiency on account of the somewhat decreased rate of 
discharge, equation (3) remains unchanged. Equation (4) 
passes over into 

Pa = Po + ~ V <P- a " 1) O) 

to which must be added the continuity equation 

P 0 v x =• P 2 v 2 (10) 

and finally the energy equation in the following form. 

The difference in the total energies at two cross 
sections immediately in front of and behind the radiator 
block is equal to the quantity of heat transmitted by the 
radiator to the air. The frontal and frictio'nal resist- 
ance of the radiator gives rise to no loss in the total 
energy since the work required is directly converted into 
the friction lo'ss: 

P Q V, J gf(i 2 + A ~ :H + A ^)] = Q. 

For the heat content, we have in general 



(11) 



2 Continued from page 7 

from that of the cold radiator. A similar consideration 
of the momentum was given by H. Lorenz for the ducted 
radiator and appreciable deviations were obtained in the 
results of the measurements, ascribed by Lorenz to the 
simplified assumptions made. This contradiction is rather 
to be explained by the fact that through an oversight the 
increase in momentum was set equal to -g(P a v 2 2 - P 0 v i 3 ) an( l 
not to P 2 v s 2 - PqV! 2 which is the proper expression . (ref- 
erence 3) . 
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i -s , — a -i- ( A - kcal/mkg A' 



The quantity of heat Q(kcal/s) may be taken as a frac- 
tion of the engine power: 

Q = A a 75 N (for water cooling a ~ 0.50) 

Equation (11) then "becomes 

p 0 v J 1 —!L-(l* _ Li} + 1 ( Va a . Vx a ) ] = 75 a v ( lla ) 
L K _ 1 > p a p Q ' . 2 J 

The solution of the system of equations (3), (8a), 
(9), (10), and (lla) gives a relation between the inlet 
and outlet velocities where in place of v 1 and v 
there are again introduced the nondimens i onal velocities 

tPj = v i / v q and ^2 = v s Ao : 

n Di + cp^ (2 - n Di - c f) 



(12) 



Finally, there is obtained an equation for the velocity 
ratio cp x 

p* + 



cp 



1 L « - 1 2A 2 . ~W \ 

+ 1 + t) 

, %i .+ ^ 8 ( 2 - ^i - t ) 



+ 111 



K - 1 



+ cp s f _JL_ -n Di - l) 1! (p* + ^ ) I =6- (13) 

1 \ k - 1 1 ' k _ 1 .. . r.i. J- . ■ 1 



where 



75 a N 
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Equation (13) can "be represented in the form of a family 
of curves with the parameters , c w , j3 2 , and 5. It 

is readily shown that the two middle terms of the expres- 
sion in "brackets, which arise essentially from the squared 
velocities in equation (lla) , are' small under normal con- 
ditions compared with the other two terms. Furthermore, 
in the fourth term of the expression in brackets T)n may 

"be neglected in comparison with p* . The resulting expres- 
sion, neglecting a small term under the root, then becomes 



cp. = cp k 



1 - I 



K - 1 



j3 + 1 + 



°w 

T 75 



a N 



V 



Po F v o 



(14) 



where t\q^ 
tion (7) or 



for a given case is to be determined by 

estimated. The remaining terms on 
hand side of the equation may be considered as 
velocity through 
quantity of 
pared with 



the radiator decreases 
heat carried 'off a N . The 
the cold radiator amounts to 



Acp 



equa- 

the right- 
known. The 
linearly with the 
decrease as com- 



(15) 



The internal drag of the cold radiator is given by 
equation (l). In the same manner there is obtained for 
the internal drag of the heated radiator 



>o * v x 



(v< 



P v.) = P° F v ( 



where cp 2 is determined from equation (12). 
in drag therefore amounts to 



<P X (1 - P <Vs) 

The decrease 



= ff ik C 1 - ~) .= (i - — -; w i> ( i6 > 

\. . Cp k 1 - P ^ 



Only the internal radiator drag can be readily comput- 
ed. If the change in the rate of flow in heating the radi- 
ator is not very large, the external drag may be considered 
as unchanged. Equation (16) in that case gives, to a good 
degree of approximation, the total decrease in the drag due 
to the heating. 

In some cases a different presentation of these re- 
sults is desirable. The heat given off by a radiator ele- 
ment is generally expressed in the following form 
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Q = 3* Vj g P 0 c p -nth 6 (kcal/s) . 
or "" 

Q = 75 a N = 427 I 1 v Q g P Q c p -n th 6 q> x (mkg/s) . 

where T,^ is 'the thermodynamic efficiency of the radi- 
ator, which is known from tests on the radiator Mock in 
the duct as a function of the discharge rate, 6 is the 
temperature drop between the cooling fluid and the air at 
the radiator inlet. With this expression for Q or a IT 
the simplified equation' (13) gives 



cd 1 = 



CP-. 



J\ + c e~ 



(17) 
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ttj-i s p o g p "nth P + 1 + "2" 



j3 a + c w -(1 - 



or with 



K = 1.4, c p = 0.24 

P 0 p + 1 + 

Po u " (3 3 + c w -(1 - 7W) 



C = 288 ^ Tl th ~ (17a) 



where, except for p" and , C contains only magni- 

tudes that are known from. experiment for a ducted radiator. 
If, in addition, there is known the velocity cp k for the 
cold ducted radiator r,j may he determined from equation 
(7) and the velocity for the heated radiator condition (e) 
from equations (17) and (17a). The exit velocity cp 2 , the 
decrease in the discharge velocity and in the drag are giv- 
en by equation (12), (15), and (15), respectively. 

By the method described above, the decrease in the 
discharge rate and in the internal drag was computed for 
radiators of various drag coefficients and thermodynamic 
efficiencies for 6 = 55° (mean value for water cooling) 
and for 6 = 115°, TU being estimated at 0.85 and 
p 



P 



-2. = 
o 



0.0000135 corresponding to normal conditions at 



3,000 meters altitude. The results are presented in fig~ 
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ures 2 to 5 as a function of- the cross-section ratio F a /F 

for Ti tll = O.S and 0.. 7 . As may be seen for 6 = 115° 

ri,. = 0.7 and F Q /F = 0.15, the condition is reached at 

'tn a v 

which the radiator no longer possesses any internal drag. 
Since, however, for such duct openings the internal drag 
(amy amount to only a. very small portion of the total drag, 
the decrease based on the total drag is relatively small 
unless special means are provided for keeping the external 
drag small. The decrease in the discharge rate is prac- 
tically independent of c w and in the figures is repre- 
sented therefore by only one curve. The effect of the de- 
crease in the radiator drag on the power required for tow- 
ing the duct system may readily be determined by the usual 
method and will therefore be' omitted. 



Translation by S. Reiss, 
National Advisory Committee 
for Aeronautics. 
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Figure 1.- Flow through and around a ducted radiator. 
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Fig. 2 
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Fig. 4 
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Fig. 3 
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Fig. 5 
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Figures 2 to 5.- Decrease in the discharge velocity and in the internal 

drag of a heated radiator in comparison with the unheated 
Condition as a function of the duct opening ratio F a /F for various drag 
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